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Distance to G14.33— 0.64 in the Sagittarius Spiral Arm: 
H 2 Maser Trigonometric Parallax with VERA 
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Abstract 

We report on trigonometric parallax measurements for the Galactic star forming region G14.33— 0.64 
toward the Sagittarius spiral arm. We conducted multi-epoch phase- referencing observations of an H 2 
maser source in G14.33— 0.64 with the Japanese VLBI array VERA. We successfully detected a parallax 
of 7r = 0.893 ±0.101 mas, corresponding to a source distance of d = 1.12 ±0.13 kpc, which is less than 
half of the kinematic distance for G14.33— 0.64. Our new distance measurement demonstrates that the 
Sagittarius arm lies at a closer distance of ~ 1 kpc, instead of previously assumed ~ 2 — 3 kpc from kinematic 
distances. The previously suggested deviation of the Sagittarius arm toward the Galactic center from the 
symmetrically fitted model (Taylor & Cordes 1993) is likely due to large errors of kinematic distances at low 
galactic longitudes. G14.33— 0.64 most likely traces the near side of the Sagittarius arm. We attempted 
fitting the pitch angle of the arm with other parallax measurements along the arm, which yielded two 
possible pitch angles of i = 34. °7 ± 2.°7 and i = 11. °2 ± 10. °5. Our proper motion measurements suggest 
G14.33— 0.64 has no significant peculiar motion relative to the differential rotation of the Galaxy (assumed 
to be in a circular orbit), indicating that the source motion is in good agreement with the Galactic rotation. 

Key words: Galaxy: kinematics and dynamics — Galaxy: structure — ISM: Hll regions — ISM: 
individual (G14.33-0.64) — masers (H 2 0) 



1. Introduction 

The Milky Way is known to be a spiral galaxy, and its 
structure has been intensively studied for many decades 
(e.g., Oort, Kerr & Westerhout 1958; Dame et al. 1987, 
2001). However, there is still little agreement on the de- 
tailed spiral pattern, including the number of the spiral 
arms (e.g., Cohen ct al. 1980; Drimmel 2000; Russeil 2003; 
Benjamin et al. 2005; Dame & Thaddcus 2008; Hou et al. 
2009). Spiral arms are regions of active star formation 
and traced primarily by Hll regions, where young stellar 
populations (hot OB stars) ionize surrounding gas. The 
major difficulty in revealing the precise spiral structure of 
the Galaxy arises from the lack of accurate distances to 
the Hll regions. 

Optical distance measurements such as can be obtained 
from photometric studies are limited in the Galactic disk 
by the large opacity due to dust. Instead, the most widely 
used method to map the Galaxy is to adopt kinematic dis- 
tances, which are derived by matching the observed radial 
velocities (obtained from the Doppler shift in observed fre- 
quencies) with respect to the local standard of rest (LSR) 
with line-of-sight velocities expected from a Galactic rota- 
tion model (e.g., Schmidt 1965; Brand & Blitz 1993). The 
famous work done by Georgelin & Georgclin (1976) adopts 
this method (with the help of optical observations where 



available) to map Hll regions in the Galaxy. However, 
significant unmodcllcd deviations from circular motions 
can cause large distance errors (Burton & Bania 1974). 
Accurate and direct distance measurements without any 
assumption on the Galactic rotation are thus of the great- 
est importance to delineate the true Galactic structure. 

It has become feasible to map Galactic structure with 
VLBI (Very Long Baseline Interferometry) techniques, no- 
tably with the phase-referencing VLBI technique, by di- 
rectly measuring trigonometric parallaxes of strong maser 
sources in star- forming regions associated with Hll re- 
gions throughout the Galaxy. In addition to precise dis- 
tances and absolute sky positions that locate the source 
in 3 dimensions in the Galaxy, measurements of absolute 
proper motions yield the full 3-dimcnsional space mo- 
tions (i.e., secular proper motions and source distances 
together give tangential velocities), which enables one 
to obtain full source information for Galactic structure 
and dynamics. Reid et al. (2009b) recently refined our 
knowledge of the Galactic spiral structure and kinemat- 
ics by integrating early results from VLBI astrometry of 
the Galaxy for total 18 high-mass star-forming regions 
(HMSFRs) with methanol, H2O, SiO maser and contin- 
uum emission, carried out with the NRAO Very Long 
Baseline Array (VLBA) and with the Japanese VERA 
project. VERA (VLBI Exploration of Radio Astrometry) 
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Table 1. VERA Observations of G14.33-0.64 



Epoch 
(1) 


Date 
(2) 


Day of Year 
(3) 


Time Range [UT] 
(4) 


Beam [ 
(5) 


mas] 


Beam E l>35° 
(6) 


[mas] 


1 


2006 Oct 27 


2006/300 


03:00-12:00 


1.87x0.89 ( 


a -25. °5 


2.70x0.75 @ 


-7.°3 


2 


2006 Nov 26 


2006/330 


01:00-08:45 


1.81x0.83 <S 


a -26. °7 


2.78x0.73 @ 


-5.°9 


3 


2007 Jan 7 


2007/007 


22:00-05:45 


1.87x0.86 ( 


a -26. °i 


— 




4 


2007 Feb 14 


2007/045 


20:00-03:43 


1.92x0.86 ( 


a -24.°2 


2.48x0.81 @ 


-4.°5 


5 


2007 Mar 27 


2007/086 


17:00-00:43 


2.10x0.82 ( 


a -28. °3 


2.67x0.76 @ 


-5.°8 


6 


2007 May 6 


2007/126 


14:00-21:43 


2.24x0.82 ( 


a -26. °0 


2.49x0.82 @ 


-8.°1 


7 


2007 Aug 8 


2007/220 


08:00-15:50 


(1.82x0.92 ( 


a -22. °0) 


2.64x0.77 @ 


-l.°6 


8 


2007 Oct 10 


2007/283 


04:00-11:55 


(1.72x0.92 ( 


3 -24.°6) 


2.72x0.75 @ 


-3.°5 


9 


2008 Jan 16 


2008/016 


21:00-04:55 


(1.79x0.92 ( 


a -27.°5) 


2.49x0.81 @ 


-6.°4 


10 


2008 Apr 14 


2008/105 


15:00-22:55 


(2.00x0.85 ( 


a -30.°9) 


2.69x0.75 @ 


-8.°6 


11 


2008 Jul 21 


2008/203 


08:30-16:15 


(1.89x0.84 ( 


a -24.°5) 


3.03x0.70 @ 


-8.°1 


(1) Epoch number. (2) The date of observation start time in universal time (UT). (3) Day ot year of observation. 
(4) Start time and end time in UT. (5) Beam size (major and minor axes) and its position angle (PA) east of 
north in single-beam images (with no data flagged). Parentheses indicate epochs not used in relative proper- motion 
measurements since the reference spot 4b and feature 4 were not detected. (6) Beam size and its PA east of north 
in dual-beam phase-referenced images, where data with antenna elevations below 35° were nagged. 



is the first VLBI array dedicated to phase referencing 
VLBI for Galactic astrometry, consisting of 4 antennas 
(20 meters each in diameter) across Japan (e.g., Honma 
et al. 2000). The recent VERA results for Galactic as- 
trometry through maser parallax measurements are re- 
ported by Honma et al. (2007), Hirota et al. (2007, 2008a, 
2008b), Imai et al. (2007), Sato et al. (2008), Kim et al. 
(2008), Choi et al. (2008), Nakagawa et al. (2008) and Oh 
et al. (2009). 

The object of this study, G14.33-0.64 (IRAS 
18159—1648), is a Galactic star-forming region and is 
VERA's first target source toward the Sagittarius spiral 
arm in the inner Galaxy, which is an important step to- 
ward our goal of mapping the structure of the Galaxy. 

In particular, located at a low galactic longitude of I = 
14.°33 (with a latitude of b = -0.°64 within the Galactic 
plane), G14.33— 0.64 is expected to trace the closest part 
of the Sagittarius arm to the Sun, and thus is an important 
target to determine the direct distance to the arm. 

G14.33— 0.64 was initially discovered as a far-infrared 
(FIR) source in a 70- /^m survey of the Galactic plane by 
Jaffe, Stier, & Fazio (1982). It was soon followed by the 
first detection of H2O maser emission at 22 GHz associ- 
ated with the FIR source by Jaffe, Glisten, & Downes 
(1981). Later the H2O maser emission was identified 
with an IRAS point source by Scalise, Rodriguez, & 
Mendoza- Torres (1989). Both class I and II methanol 
(CH3OH) maser sources were also found in the region: 
class II emission at 6.7 GHz (Walsh et al. 1995, 1997) 
and class I emission at 36 GHz, at 44 GHz (Slysh et al. 
1994, 1999), at 84 GHz (Kalenskii et al. 2001), and at 
95 GHz (Val'tts et al. 2000). G14.33-0.64 has been ob- 
served to display OH thermal absorption line at 1665 
MHz (Wouterloot et al. 1993), NH 3 (1,1) and (2,2) inver- 
sion transition lines at 23.7 GHz (Molinari et al. 1996), 
CS( J = 2 — > 1) and CS(J = 5—5-4) rotational transition 
lines at 98.0 GHz (Bronfman et al. 1996) and at 244.9 GHz 
(Shirley et al. 2003), respectively, and 1.2-mm continuum 



emission (Faundez et al. 2004). The radial velocities ob- 
served for many molecular lines of G14.33— 0.64 are in 
good agreement at Vlsr — 22 km s _1 . 

In the present study, we report on our successful deter- 
mination of the parallax of G14.33— 0.64 with VERA as 
a step toward revealing the structure of the Sagittarius 
spiral arm in the inner Galaxy. 

2. VERA Observations 

VERA observations of the 22-GHz H2O maser source 
(the 616 — »523 rotational transition) in G14.33— 0.64 were 
carried out at 11 epochs between 2006 October and 2008 
July as listed in Table 1. Using VERA's dual-beam 
mode for phase referencing (e.g., Honma et al. 2003, 
2008a), we simultaneously observed the H2O maser source 
in G14.33— 0.64 and the extragalactic position-reference 
quasar (phase calibrator) J1825— 1718 with an angular 
separation of 1°.7 at a position angle (PA) of 108° east 
of north relative to G14.33— 0.64. The flux density of the 
phase calibrator J1825— 1718 was up to ^140 mJy. A 
nominal H2O maser position for G14.33— 0.64 was used 
as reference center both for observation and for cor- 
relation: a 20 oo =18 h 18 m 53. s 8 and <5 2 ooo = -16°47'50"0 
(Comoretto et al. 1990). The position of J1825-1718 
was adopted from the second VLBA Calibrator Survey 
by Fomalont et al. (2003): a 20 oo =18 h 25 m 36. s 532283 and 
£2000 = -17°18'49."84781. The ICRF source NRAO 530 
(J1733— 1304; Ma et al. 1998) was also observed as a bright 
calibrator source (fringe finder) for 7-minutc scans hourly 
in each beam. 

The instrumental phase difference between the two 
beams were calibrated by recording the real-time phase 
data with an artificial noise source in each beam 
(Kawaguchi et al. 2000; Honma et al. 2008a). Left-hand 
circularly polarized signals were digitized at 2-bit sam- 
pling and recorded at a data rate of 1024 Mbps. In the 
total bandwidth of 256 MHz (16x16 MHz), one of the 
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sixteen 16-MHz IF channels was assigned to the H2O 
maser lines in G14.33— 0.64. The other 15 IF channels 
were for the continuum emission in the phase calibrator 
J1825— 1718, with the central IF channel set at the maser 
frequency, using the VERA digital filter unit (Iguchi et al. 
2005). 

The data correlation was performed with the Mitaka 
FX correlator (Chikada et al. 1991). In order to obtain 
sufficient resolution for the H2O maser lines, only the cen- 
tral 8 MHz (of the 16-MHz IF channel) for the maser lines 
was split into 512 spectral points, yielding frequency and 
velocity resolutions of 15.625 kHz and 0.21 km s _1 , re- 
spectively Due to the spectral splitting method, one of 
the other 15 IF channels for J1825— 1718 was also split 
into 512 spectral points (with the maser channel), which 
was not used for data reduction. The other 14 IF channels 
were split into 64 spectral points each and used in data 
reduction. 

The system noise temperatures at the zenith were typi- 
cally T sys =150—300 K for the first 5 epochs. For the last 
6 epochs, one or two antennas showed higher system noise 
temperatures of T sys = 300 — 800 K due to bad weather, 



i-sys 

while the other antennas remained at T< 



sys 



=150-300 K. 



Data Reduction 



Visibility calibration and imaging were performed in 
a standard manner with the NRAO Astronomical Image 
Processing System (AIPS) package (Greisen 2003). The 
observed frequencies of the maser lines were converted to 
radial (line-of-sight) velocities with respect to the local 
standard of rest (LSR), Vlsr, using a rest frequency of 
22.235080 GHz (Pickett ct al. 1998) for the H 2 6i 6 ^5 23 
transition. 

We first searched for the relative positions of all H2O 
maser spots in the single-beam data (i.e., without phase- 
referencing to the calibrator J1825— 1718 in the other 
beam) of the third epoch and found maser emission 
over several spectral components (see figure 1). At this 
epoch, the brightest H 2 maser channel was at Vlsr = 
26.6 km s _1 (feature 7 in table 3), and the visibilities 
of all maser channels were firstly phase-referenced to this 
channel by fringe fitting (AIPS task FRING) using the 
channel and by applying the phase solutions to all the 
maser channels. In order to find the maser spot distribu- 
tion, we imaged all channels with the AIPS task IMAGR 
with a wide field of view of ~ 2" x 2" around the reference 
maser spot (feature 7), with 2048 pixels x 2048 pixels of 
size 1 mas. Many of the maser spots were outside this 
field (~ 5" offset from feature 7 as seen in table 3) and 
were found by fringe rate mapping with the AIPS task 
FRMAP and by shifting the image center accordingly. 

3.1. Phase Referencing for Parallax Measurements 

Next, we obtained absolute-position maps of 
bright maser spots by phase referencing visibilities 
of G14.33-0.64 to those of J1825-1718. For each epoch, 
phase solutions from fringe fitting with J1825— 1718 were 
applied to the H 2 maser channels of G14.33— 0.64 for 
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Fig. 1. Spectral evolution of H2O maser emission in 
G14.33— 0.64. Numbers show the observed year and day 
of year. Scalar-averaged cross-power spectra are shown be- 
tween VERA Mizusawa and Iriki stations. The radial ve- 
locities for many molecular lines of G14.33— 0.64 arc at 
l/ LSR ~22 km s- 1 . 

the corresponding frequencies. The instrumental phase 
difference between the two beams was also corrected using 
the real-time phase-calibration data recorded during each 
observation. Visibility phase errors caused by the Earth's 
atmosphere were calibrated based on GPS measurements 
of the atmospheric zenith delay which occurs due to 
tropospheric water vapor (Honma et al. 2008b). 

Since a nominal reference center of G14.33— 0.64 was 
used for correlation, we first imaged the phase-referenced 
maser data to find the positional offset of each maser 'fea- 
ture' (i.e., a group of maser spots in the same position 
over adjacent velocity channels) from the reference center, 
and then recalculated and corrected the delays for the ob- 
tained absolute positions of the maser features until the 
features came at the map center within 10 mas. After 
correcting the absolute position of the reference center for 
each maser feature, we used the same map center at all 
epochs for the same maser feature. We imaged the de- 
tected maser spots with the AIPS task IMAGR for a field 
of view 25.6 mas x 25.6 mas around each map center, with 
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512 pixels x 512 pixels of size 0.05 mas. Maser positions 
were fitted with elliptical Gaussian distributions with the 
task JMFIT. RMS noise levels in each image per channel 
were typically 200—600 mJy/beam. 

We performed least-squares fitting to simultaneously 
solve for the sinusoidal parallax curve and linear proper 
motion in right ascension (RA) for maser spots at consecu- 
tive velocity channels for two features that were persistent 
over more than a year. We did not solve for the parallax in 
declination because positional errors due to tropospheric 
zenith delay residuals were larger in declination as in other 
measurements (e.g., Sato et al. 2007, 2008) and also be- 
cause angular resolution was lower in declination than in 
RA for G14.33— 0.64 (see table 1 for beam size), and the 
parallax ellipse was smaller in declination. Instead, we 
removed the parallax obtained from RA fits to fit linear 
proper motion in declination. 

Since image distortion and positional errors due to tro- 
pospheric zenith delay residuals are severe for sources at 
low elevation angles associated with low source declina- 
tions including G14.33-0.64 (e.g., Honma et al. 2008b), 
we attempted 4 different elevation cutoff values 25°, 30°, 
35°, 39°, below which we flagged the data with the AIPS 
task UVFLG. For cutoffs above 39°, imaging became dif- 
ficult with high sidclobes due to insufficient data. We 
adopted an elevation cutoff of 35° to obtain the best fit- 
ting result. For example, a typical error in the position 
measurement with one maser spot reduced from 0.18 mas 
to 0.14 mas by changing the elevation cutoff from 30° to 
35°. Flagging low-elevation data changed the beam size 
of antennas to be elongated in declination, however the 
beam size in RA was kept almost unchanged or slightly 
better (smaller) (see table 1). 

3.2. Single-Beam Analysis for Relative Proper Motions 

We also measured relative proper motions from the 
single-beam data to study internal motions of H2O maser 
spots. Since the H2O maser emission in G14.33— 0.64 was 
variable over the observing period (figure 1), the phase- 
reference maser channel used for fringe fitting differed 
epoch to epoch: we used the brightest velocity channel 
at each epoch as the phase reference, excluding the chan- 
nels around Vlsr ~ 26 km s" 1 (feature 7 in table 3) be- 
cause the maser spots in this velocity range were 5" away 
from the other spots. We imaged each maser spot with 
the AIPS task IMAGR for a field of view of 25.6 mas x 
25.6 mas (512 pixels x 512 pixels of size 0.05 mas) by 
shifting the map center. The FWHM beam size of each 
epoch is shown in table 1. RMS noise levels in each image 
per channel were typically 50—110 mJy/beam. 

Maser positions were fitted with elliptical Gaussian dis- 
tributions with the task JMFIT and were measured rel- 
ative to the reference spot chosen at each epoch. In or- 
der to obtain relative proper motions of all spots, we cal- 
culated all maser positions relative to the maser spot at 
Vlsr ~ 21.6 km s _1 (spot 4b in table 3) by subtracting 
the position of this spot from the maser positions at each 
epoch. Since feature 4 (including spot 4b) was only persis- 
tent over the first 6 epochs, relative proper motions were 



measured over the first 6 epochs. 

Our criteria for detection of a maser feature are: (1) a 
signal-to-noise ratio higher than 7 is obtained in the map 
at more than two consecutive velocity channels, (2) the 
spots are identified at three or more epochs for detecting 
relative proper motions, and (3) their positions agree with 
those expected from the fitted proper motions within 1 
mas. In table 3, we also list the strong feature at Vlsr ~ 
26 km s , even though it has no measured proper motion. 

4. Results 

Fig 1 shows the spectral evolution of H2O maser emis- 
sion in G14.33— 0.64 over the observing period: scalar- 
averaged spectra are shown for the baseline between 
VERA Mizusawa and Iriki stations. 

Parallax Measurements 

Table 2 summarizes the results from measurements of 
parallax 7r in RA (X) and absolute proper motions \xx 
and fiy in RA (X) and Dec (Y). We used a total 
of seven maser spots of two maser features (features 1 
and 4; feature IDs in table 2 correspond to those in ta- 
ble 3). The absolute maser positions used for the mea- 
surements were: 012000 =18 h 18 m 54. s 67444 and <5 2 ooo = 
-16°47'50"2640 for feature 1; a 20 oo =18 h 18 m 54. 8 65341 
and (5 2 ooo = -16°47'50."0650 for feature 4. 

Errors in the measurements are indicated in parenthe- 
ses in table 2. For single-spot measurements, errors were 
estimated from the residuals from the least-squares fit- 
ting with uniform weights for all epochs. For combined 
fits where different spots or features were simultaneously 
fitted with a single parallax and with different proper mo- 
tions, we have estimated the upper limit of the errors. We 
will discuss error estimates further in detail in §5.1. 

The final value of the parallax (from the combined fit 
with all seven spots) is tt = 0.893 ±0.101 mas. This 
corresponds to a source distance of d = 1.12 ±0.13 kpc. 
Absolute proper motions \ix and \iy of the seven spots 
listed in table 2 were derived using this final value of tt, 
instead of using different tt values from single-spot mea- 
surements. 

Figures 2 and 3 show the position measurements of the 
seven spots used for parallax and absolute proper motion 
fitting. Numbers indicate feature IDs corresponding to 
those in table 2. Figures 2a and 2b show eastward (X) 
and northward (Y) positional offsets versus time, respec- 
tively, for seven maser spots. Additional constant offsets 
are added to each maser spot in the figures for clarity. The 
best-fit models for the single parallax (solid curves) and 
different proper motions (gray lines) are plotted for the 
seven spots from the combined fit. Error bars are plot- 
ted for the standard deviation (a) of the post-fit residuals 
from the least-squares fitting. Figure 3 shows the trajec- 
tory on the sky. 

4.2. Proper Motions 

Table 3 lists the results from relative position and 
proper-motion measurements from the first 6 epochs. A 
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Table 2. Parallax fitting for G14.33— 0.64 with elevation cutoff 35 degrees. 



Feature ID 


Vlsr 


-^epochs 


Detected Epochs 


RA Parallax, ir 


P-x 


My 


# 


[km s -1 ] 






[mas] 


[mas yr _1 ] 


[mas yr _1 ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


la 


14.6 


6 


45678-10- 


0.931 (0.124) 


6.13 (0.27) 


-4.50 (0.37) 


11) 


14.8 


8 


4567891011 


0.936 (0.151) 


6.47 (0.19) 


-4.15 (0.26) 


lc 


15.0 


8 


4567891011 


0.950 (0.141) 


6.49 (0.19) 


-4.12 (0.26) 


Id 


15.2 


6 


45-78910- 


1.004 (0.135) 


6.28 (0.26) 


-4.23 (0.35) 


1 combined 








0.954 (0.130) 






4a 


21.4 


6 


12-456 10- 


0.629 (0.171) 


-1.60 (0.23) 


-0.26 (0.30) 


4b 


21.6 


6 


12-456 10- 


0.631 (0.162) 


-1.58 (0.23) 


-0.42 (0.30) 


4c 


21.8 


6 


12-45 8-10- 


0.900 (0.151) 


-1.70 (0.21) 


-0.12 (0.28) 


4 combined 








0.768 (0.160) 






1&4 combined 








0.893 (0.101) 







(1) Feature/spot ID, corresponding to table 3. (2) LSR velocity of the mascr spot. (3) Total number of detected 
epochs. (4) Detected epochs. (5) Measured parallax in right ascension in mas (with estimated errors in parentheses). 
(6) (7) Proper motions in right ascention and in declination, respectively. The results presented here were obtained 
by fitting with a single parallax of 0.893 mas (the final result from the combined RA parallax fit). 
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Fig. 3. Trajectory of maser positions on the sky. Reference 
positions are the same as in figure 2 for each feature. 
Additional shifts are given for clarity: Ay = +9, 5, 1, —3 
mas for la, lb, lc, Id; —7, —9.5, —12 mas for 4a, 4b, 4c, 
respectively. 

total of 6 maser features are presented here with proper 
motions detected over 3 or more epochs. 

Errors of relative proper motions (shown in parentheses 
for fx x and \i y in table 3) for each spot were estimated from 
formal fitting uncertainties scaled by the RMS residuals 
of the spot positions. For each feature (#1 through 6), 



the relative position and proper motion were calculated 
as error- weighted means of those for all detected spots of 
the feature, as notated by 'w-mean' in table 3. 

Figure 4 shows the maser distribution and proper mo- 
tion for G14.33— 0.64. Figure 4a and 4b are radio maps of 
the region with contours showing continuum emission at 
6-cm wavelength (C-band) from VLA archive data (pro- 
gram AH361) observed in the 'C configuration at an an- 
gular resolution of 3" (Hughes & MacLeod 1994). Our 
VLBI absolute positions of H2O maser features are also 
shown. Our absolute position accuracy is essentially lim- 
ited by the position errors of the position-reference quasar 
J1825-1718, which are 1.23 mas in RA and 1.97 mas in 
Dec (Fomalont et al. 2003). 

Hughes & MacLeod (1994) originally associated the 
H2O maser emission in G14.33— 0.64 with brighter radio 
continuum emission, offset 1.'16 toward northeast from the 
IRAS position (at the origin) as seen in figure 4a. Our 
new VLBI map finds the H2O maser emission (feature 7 
in particular) associated (within 5") with a different and 
closer radio continuum source and yields the first precise 
distribution of H2O maser spots in G14.33— 0.64. 

Figure 4c shows the absolute proper motions of maser 
features 1 to 6, which were obtained by adding relative 
proper motions (in table 3) to the absolute proper motion 
of the reference spot 4b (in table 2). These absolute proper 
motions are not corrected for apparent motions due to 
the Solar motion and the Galactic rotation, in addition 
to the peculiar motion of the source. The map origin is 
the position of the reference spot 4b at the first epoch: 
a 2 ooo =18 h 18 m 54. s 653181 and 62000 = -16°47'50."07668. 

Figure 4d shows the internal motions of the maser fea- 
tures relative to the mean motion of the features. The 
mean motion of all features 1 to 6 was obtained by aver- 
aging obtained proper motions over the 3 distinct regions: 
(1) feature 1; (2) features 2 and 3; and (3) features 4, 
5 and 6. We took an unweighted mean of relative mo- 
tions of maser features in each region, and then took an 
unweighted mean of the 3 regions, as listed in table 3 
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by '1,23,456 u-mean'. We obtained mean relative mo- 
tion of (fix, /Jj,)=(2.53, —2.08) mas yr _1 (the bar symbols 
indicate mean values). By adding the absolute proper 
motion of the reference spot 4b, (fix, HY)ib =(—1.58, 
—0.42) mas yr _1 , we obtained the absolute mean motion 
{fix, ^y)=(0.95, —2.50) mas yr -1 . Note a proper mo- 
tion of 1.00 mas yr -1 corresponds to a linear velocity of 
5.31 km s _1 at a source distance of 1.12 kpc. In §5.5, 
we will adopt this mean motion to discuss the systemic 
motion of G14.33— 0.64, by taking errors into account to 
allow for the possibility that the mean maser motion docs 
not trace the systemic motion. 

5. Discussion 

5.1. Astrometric Error Sources 

In this section, we will discuss possible error sources in 
our parallax and proper-motion measurements and how 
we estimated the errors. 

The first possible error source in individual position 
measurements is thermal errors due to noise, which can be 
approximated by the halfwidth (HWHM) of the beam size 
divided by the signal-to-noise ratio of the maser map. We 
find thermal errors can account for the errors of relative 
position measurements in the single-beam data. Thermal 
errors are ~ 0.01 — 0.1 mas (i.e., beam size ~ 1 mas, sig- 
nal ~ 1 — 10 Jy/beam, and noise ^O.lJy/beam) and agree 
well with errors in the relative proper motions as listed in 
table 3, which were estimated from standard deviations 
from the post-fit residual from the least-squares fits (see 
§3.2 and §4.2). 

However, for parallax and proper-motion measurements 
in the dual-beam data, errors in the measurements are 
larger than expected from thermal errors of ~ 0.1 mas 
(i.e., beam size ~ 1 mas, signal ^ 3 Jy/bcam, and noise 
^0.3 Jy/beam). The standard deviations from the fits 
were a = 0.26 mas and thus are larger than thermal noise 
errors. 

Here we do not consider the reference quasar as pre- 
dominant error source since it did not show any resolved 
structure. Also, even though the accuracy of the maser 
absolute position is limited by the uncertainties of the 
reference quasar position, the positional error of the ref- 
erence quasar only adds as a constant offset to the maser 
spot position at each epoch and do not contribute to un- 
certainties in parallax and proper-motion measurements. 

One of the likely sources that would cause large errors 
in the parallax and proper-motion measurements is errors 
in modeling of troposphcric zenith delay (see Sato ct al. 
2008 and references therein). Indeed, the fact that a high 
elevation cutoff of 35° yielded the best-fit result for the 
parallax fitting for G14.33— 0.34 indicates that this low- 
declination source is subject to troposphcric delay errors. 
However, if errors in modeling of troposphcric zenith delay 
are the predominant error source, then all maser features 
at the same epoch should show systematic errors in the 
position measurements. As can be seen in figure 2, the 
deviations from the parallax fits clearly differs for the two 
different maser features at each eposh (features 1 and 4), 



which indicate that errors are random for different fea- 
tures at the same epoch. Therefore, it is likely that errors 
in modeling of troposphcric zenith delay are not the pre- 
dominant error source in the remaining data after having 
removed as much effects of the tropospheric delay errors 
as possible by adopting a high elevation cutoff. 

Another likely error source in the parallax measure- 
ments is variation in maser structure. In figure 2, the 
tendency of deviations from the parallax fits is similar for 
maser spots in the same feature but different between dif- 
ferent features (features 1 and 4). This is consistent with 
the fact that the variation of maser structure causes po- 
sitional errors that are uncorrelated for different features 
but might be correlated between maser spots in adjacent 
velocity channels within the same feature. In our parallax 
measurements, the variation of maser structure is likely 
the predominant error source. 

We estimated errors of the parallax measurements from 
the post-fit residuals from the least-squares fitting. For 
different spots within the same maser feature (e.g., spots 
la, lb, lc, and Id in feature #1), we allow for the pos- 
sibility that errors due to variation in maser structure 
data may be partially correlated. As a conservative ap- 
proach, wc assumed errors of all spots within the same 
maser feature at the same epoch are 100% correlated (but 
random for different features). This means that, even 
though we used 7 maser spots of la, lb, lc, Id, 4a, 4b, 4c 
for the measurements, we assume only 2 different maser 
features contribute as 2 independent spots to reduce the 
errors of the measurements. We obtained an error of 
ov = 0.101 mas. Instead, if we assume errors due to vari- 
ation in maser structure are random and uncorrelated for 
all the 7 spots (la— Id and 4a— 4c), the errors would reduce 
to a'^ = 0.060 mas. In reality, errors of different spots in 
the same feature are not likely to be 100% correlated but 
only partially correlated (if not uncorrelated). Therefore, 
the error estimate of ov =0.101 mas in our parallax mea- 
surements is the upper limit of the errors, adopted as a 
conservative approach. 

5.2. Distance to the Sagittarius Spiral Arm 

Our parallax measurement for G14.33— 0.64 reveals the 
source distance to be d= 1.12±0.13 kpc, which is less than 
half of previously derived kinematic distances. The kine- 
matic distances for G14.33— 0.64 are, for example, 2.5 kpc 
by Molinari et al. (1996) from the NH 3 (1,1) and (2,2) 
lines; 2.6 kpc both by Walsh ct al. (1997) and Val'tts et al. 
(2000) from 6.7-GHz and 95-GHz methanol maser lines, 
respectively; 3.1 kpc from the HllOa line and 2.6 kpc 
from H2CO absorption lines by Sewilo et al. (2004). All 
of the kinematic distances above were derived using the 
Galactic rotation model by Brand & Blitz (1993). Palagi 
et al. (1993) derived 2.7 kpc from H2O maser lines with 
a peak at Vlsr = 22.8 km s _1 using the rotation curve 
of Brand (1986). The good agreement among previous 
kinematic distances is a result of using the same rotation 
model and similar radial velocities Vlsr — 22 km s _1 ob- 
served at different wavelengths. The most persistent H2O 
maser feature in our measurements, feature 4, also showed 
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Fig. 7. Galactic maser source locations in the Sagittarius 
and Local (Orion) arms, superimposed on artist's con- 
ception (R. Hurt: NASA/JPL-Caltech/SSC). '0' indi- 
cates the location of the Sun and '*' the position of the 
Galactic center. The red square shows the new loca- 
tion of G14.33— 0.64 based on our parallax measurements. 
Three star-forming regions, G35.20— 0.74 (blue diamond), 
G35.20-l.74 (pink triangle) and W51 IRS2 (green hexagon), 
possibly belonging to the Sagittarius spiral arm, are also 
indicated with parallactic distances. The positions of 
five sources in the Local (Orion) 'arm' or spur are indi- 
cated by cyan dots with precise parallactic distances (see 
text). Errors for all parallactic distances are also shown, 
which arc mostly smaller than the size of the symbols. 

a radial velocity of Vlsr — 22 km s _1 , which agrees well 
with the systemic radial velocity of G14.33— 0.64, but sev- 
eral other spectral components differed up to 10 km s _1 
in radial velocities. 

Figure 5 shows the classic model of the Galaxy by 
Gcorgelin & Georgelin (1976). Gray lines show the 
modified model by Taylor & Cordes (1993). Note that 
a shift toward the Galactic center in the position of 
the Sagittarius arm was introduced by Taylor & Cordes 
(1993), to correspond better with the kinematic distances 
of Downes et al. (1980). Dowries et al. (1980) estimated 
the kinematic distances to Galactic H II regions from radio 
observations of HI 10a and H2CO lines using the Schmidt 
(1965) model, with typical errors of ±1 to 2 kpc for galac- 
tic longitudes I = 20° to 60°, which can be more than 
~50% errors for the Sagittarius arm at lower galactic lat- 
itudes. Although G14.33— 0.64 was not in the catalog by 
Downes et al. (1980), it was in the catalog by Sewilo et al. 
(2004) in HI 10a and H2CO line observations. It can be 
clearly seen in figure 5 that the kinematic distance (shown 
as the yellow square) places G14.33— 0.64 as well as the 
interpolated Sagittarius arm further toward the Galactic 
center, like other sources in the arm. 

However, our direct parallax measurements (red square 



in figure 5) reveals the location of G14.33— 0.64 to be 
closer to the Sun and outward the Galaxy, in good agree- 
ment with the Sagittarius arm originally modeled by 
Gcorgelin & Georgelin (1976), without the 'bump' toward 
the Galactic center. In figure 5, three other star-forming 
regions, G35.20-0.74 (blue diamond), G35.20-l.74 (pink 
triangle) and W51 IRS2 (green hexagon), possibly in the 
Sagittarius spiral arm, are also plotted with parallax dis- 
tances of 2.19±°;|^ kpc, 3.27±^;21 kpc (Zhang et al. 2009) 
and 5.1±?;| kpc (Xu et al. 2009) from the VLBA for 12- 
GHz methanol maser emission. 

It is most likely that the 'bump' in the Sagittarius spi- 
ral arm toward the Galactic center suggested in Taylor 
& Cordes (1993) is due to errors of kinematic distances. 
A more recent model by Cordes & Lazio (2002), which is 
built upon the Taylor & Cordes (1993) model, also retains 
the 'bump' of the Sagittarius arm toward the Galactic 
center. Both Taylor & Cordes (1993) and Cordes & Lazio 
(2002) give models for the distribution of free electrons in 
the Galaxy, upon which most pulsar distances are deter- 
mined using the observed dispersion measures (DM), i.e., 
the column density of electrons toward the pulsars (Frail 
& Weisberg 1990). These models are built by numeri- 
cally fitting predicted and observed dispersion measures 
for pulsars with known 'independent distance estimates' 
(Taylor & Cordes 1993), most of which come from uncer- 
tain kinematic distances. 

In particular, kinematic distances are more severely af- 
fected by errors of the radial velocities for sources at low 
galactic longitudes than at high longitudes. For exam- 
ple, for the simplest assumption of circular Galactic rota- 
tion with a source distance d in the solar neighborhood 
(d <C Ro, where Rq is the distance to the Galactic center 
from the Sun), the kinematic distance d can be approxi- 
mated by c?kin ~ Vlsr/ (Asm(2l)) using Oort's constant A 
(see e.g., Karttunen et al. 2007). Errors in the kinematic 
distances (Td kin arc thus proportional to the errors in the 
radial velocities divided by sin(2^): er<i kin oc ov LSR /sin(2?). 
Therefore the kinematic distances toward the Sagittarius 
arm in the inner Galaxy are expected to be particularly 
uncertain. 

Taylor & Cordes (1993) acknowledge that pulsar dis- 
tances derived from previous models generally tend to 
be overestimated for |^| < 30° and underestimated for 
I = 50° — 70° (although they claim their own model has 
no significant dependence of distance errors on /), which 
can account for the 'bump' of the Sagittarius arm toward 
the Galactic center at low galactic longitudes. Our results 
as shown in figure 6 indicate that the previously expected 
'bump' in the Sagittarius arm toward the Galactic cen- 
ter is most likely due to errors that arise from kinematic 
distances. 

Russeil (2003) points out that the nearest part of the 
Sagittarius arm is placed at ~ 2 kpc based on kinematic 
distances (using the rotation curve of Brand & Blitz 1993), 
while a fitted regular logarithmic arm, also based on kine- 
matic distances, passes at ~ 1 kpc, indicating the pos- 
sibility that the Galaxy does not have a regular design. 
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However, our parallax measurements suggest that the 
nearest part of the Sagittarius arm, indeed, lies at ~ 1 kpc. 
The disagreement between the arm fitting and the kine- 
matic distance is likely due to errors of kinematic dis- 
tances, rather than an irregular design of the Sagittarius 
arm. 

Direct determination of distances are of great impor- 
tance and required to obtain a true map of the Galaxy 
and, in particular, of the Sagittarius arm. Our parallax 
measurement of G14.33— 0.64 with VERA reveals the lo- 
cation of the Sagittarius arm to be closer to the Sun than 
previously thought. 

5. 3. Pitch angle of the Sagittarius arm 

We attempted to fit the pitch angle i of the Sagittarius 
arm using our parallax measurement of G14.33— 0.64 with 
three other parallax measurements of sources shown in fig- 
ure 5, which may lie in the Sagittarius arm: G35.20— 0.74, 
G35.20-l.74 (Zhang ct al. 2009) and W51 IRS2 (Xu et al. 
2009). The pitch angle i is defined as the angle between 
the arm and the tangent to a Galactocentric circular orbit. 
For an ideal logarithmic spiral arm, it can be expressed 
as, ln(i?i/i?2) = — (Pi — /?2)tani, for two sources 1 and 
2 (indicated by subscripts) in the arm, where R is the 
Galactocentric radius at a Galactocentric longitude (3 (0 
toward the Sun and increasing with galactic longitude; see 
Reid et al. 2009b). 

Figure 6a shows a plot of log 10 (i?/kpc) vs. (3 (in de- 
grees) for G14.33-0.64 (red square), G35.20-0.74 (blue 
diamond), G35.20-l.74 (pink triangle) and W51 IRS2 
(green hexagon). Here we adopted the Sun-center distance 
of Rq = 8.5 kpc. Errors are indicated for each source with 
parallax uncertainties of ±lcr from this study, Zhang et al. 
(2009) and Xu ct al. (2009) . We attempted linear least- 
squares fitting to the sources with unweighted straight 
lines. (Note that we need to express lni? in natural loga- 
rithm and (3 in radians to calculate the pitch angle.) 

As seen in figure 6a, the four sources do not lie in a 
straight line, and we attempted fitting with two possible 
combinations of three sources, which are shown in gray 
lines in the figure. Line A shows a best-fit straight line for 
G14.33-0.64, G35.20-0.74 and G35.20-l.74 (excluding 
W51 IRS2), which yields a pitch angle of i = 34.°7±2.°7. 
Line B is a fitting result from G14.33-0.64, G35.20-0.74 
and W51 IRS2 (excluding G35.20-l.74), which yields a 
smaller pitch angle of i = 11. °2 ± 10. °5. This pitch angle 
i ~ 11° agrees well with the four- arm Milky Way model by 
Vallee (1995) with a best-fit pitch angle of i = 12.° 1 ± 1. 

Figure 7 is a plot of the positions of the four sources 
superimposed on artist's conception of the Milky Way. 
For comparison, five sources in the Local (Orion) 'arm' 
or spur are also shown with precise parallax measure- 
ments: G59.7+0.1 (Xu et al. 2009), Cep A (Moscadelli 
et al. 2009), Orion (Hirota et al. 2007; Menten et al. 2007; 
Kim et al. 2008), G232.6+1.0 (Reid ct al. 2009a), and 
VY CMa (Choi et al. 2009; Reid et al. 2009c). With the 
five sources, Reid et al. (2009b) fitted the pitch angle of 
the Local arm to be 27. °8 ± 4.°7, which is larger than 
the pitch angles they fitted for another spiral arm, e.g., 



16.°5±3.°1 for the Perseus spiral arm. 

In figure 6b, we also attempted a straight line fitting 
(line C) to the five Local arm sources (marked by cyan 
dots) plus W51 IRS2 (green hexagon), which yields a pitch 
angle of 26. °1 ± 12. °3, which is consistent with the pitch 
angle fitted with only five sources above. Thus the Local 
arm/spur may branch from the Sagittarius arm near the 
position of W51 IRS2, which is often considered to be at 
the tangent point of the Sagittarius arm. One possible 
interpretation is that the Sagittarius arm bifurcates near 
the position of W51 IRS2 into the Local spur (line C) at 
a pitch angle of i ~ 26° and into the other arm traced by 
G14.33-0.64 and G35.20-0.74 (line B) at a pitch angle 
of i ~ 11°. Another possibility is that the Sagittarius arm 
is traced by G14.33-0.64, G35.20-0.74 and G35.20-l.74 
(line B) and branches from the interior (Scutum-Crux) 
arm at a large pitch angle of i ~ 34°. However, more 
sources with precise parallaxes are needed to establish 
clear spiral arm structure. Ongoing and future parallax 
measurements with VERA and with the VLBA are ex- 
pected to reveal the structure of the Sagittarius arm and 
other spiral arms of the Galaxy further in detail. 

5.4- Magnetic Field Reversals and the Sagittarius Arm 

It is of interest to compare our results for the distance to 
the Sagittarius arm with studies of Galactic magnetic field 
reversals. The Galactic magnetic field has been probed 
most often by Faraday rotation measure (RM) observa- 
tions of linearly polarized emission from both pulsars (e.g., 
Noutsos et al. 2008) and extragalactic radio sources (e.g., 
Brown ct al. 2007). A common conclusion in many pul- 
sar polarization studies is that the magnetic field in the 
Local arm is clockwise while it is counterclockwise in the 
first quadrant (0° < I < 90° ) component of the Sagittarius 
arm, indicating the existence of a magnetic field rever- 
sal between the arms. Weisberg et al. (2003) found from 
pulsar polarimetry a null in the magnetic field of a width 
less than 0.5 kpc extending from near the Sun over 7 kpc 
toward / ~ 60° (figure 4 in Weisberg et al. 2003), located 
midway between the Local and Sagittarius arms, which is 
most likely the field reversal region. 

Weisberg et al. (2003) noted a "1-kpc wide strip" of 
steady magnetic field from the local reversal (midway be- 
tween the Local and Sagittarius arms) into the Sagittarius 
arm, based on the Sagittarius arm model by Cordes 
& Lazio (2002). As previously discussed, our parallax 
measurements demonstrate the Sagittarius arm lies at a 
closer distance of <~ 1 kpc, instead of previously estimated 
~ 2 — 3 kpc from kinematic distances, and we find that 
G14.33-0.64 (this study) and G35.20-0.74 (Zhang et al. 
2009) trace out the near side of the Sagittarius arm, which 
lie outside of the 'bump' delineated in Taylor & Cordes 
(1993) as well as in Cordes &: Lazio (2002). Our parallax 
measurements thus indicate that the strip of steady mag- 
netic field found by Weisberg ct al. (2003) is likely in the 
Sagittarius arm, rather than in an inter-arm region exte- 
rior to the arm. This lends support to the fact that the 
magnetic field in the Sagittarius arm is steadily and dom- 
inantly counterclockwise, and is further evidence for the 
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conclusion of Weisberg et al. (2003) that the field maxima 
tend to lie along the spiral arms, while the field reversals 
occur between the arms. 

5.5. Motion of Gl^. 33— 0.64 an d ^ e Galactic Rotation 

As seen in figures 4c and 4d, the internal motions of 
the H2O masers in G14.33— 0.64 show a bipolar jet-like 
motion on the sky, with deviations of ~ 1 — 2 mas yr -1 
from the mean, which correspond to a linear velocity of 
5 — 10 km s _1 at a distance of 1.12 kpc. The central 
radial velocity Vlsr — 22 km s _1 of the maser emission 
agrees well with other molecular line velocities, and the 
deviations up to 10 km s _1 from the central radial velocity 
agree with the proper motions. 

From the parallax, proper motion, radial velocity and 
the sky position of the H2O maser source, we can now 
calculate the full three-dimensional position and velocity 
of the source in the Galaxy. By following the methods 
described in detail by Reid et al. (2009b) to convert from 
the heliocentric reference frame to a reference frame that 
rotates with the Galaxy, we obtain the peculiar motion of 
the source with respect to the Galactic rotation. 

We adopt the mean absolute proper motion (the ref- 
erence frame in figure 4d) of (f£~x, /iV) = (0.95, —2.50) 
mas yr _1 as the systemic motion of the source (before 
the correction of the solar motion and the Galactic rota- 
tion), with uncertainties of ±2 mas yr _1 ~ ±10 km s" 1 in 
each of the eastward (X) and northward (Y) directions to 
allow for the possibility that the mean maser motion docs 
not trace the systemic motion. 

For the radial velocity, we adopted Vlsr = 22 ± 
10 km s _1 . Adopting the Hipparcos solar motion val- 
ues of Uq — 10.0 ± 0.36 km s _1 (radially toward the 
Galactic center), V@ = 5.25 ±0.62 (in the local direction 
of Galactic rotation) and Wq = 7.17 ±0.38 km s _1 (ver- 
tically upwards, i.e., toward the north Galactic pole per- 
pendicularly to the Galactic plane) from Dehnen & Binney 
(1998) with the recent best- fit results for the Galactic con- 
stants of Ro = 8.4 ±0.6 kpc and 8 = 254 ± 16 km s" 1 
by Reid et al. (2009b), and assuming a flat rotation of 
the Galaxy (i.e., rotational velocity 8 at the source lo- 
cation is the same as at the Sun, 8 ~ 80) the peculiar 
velocity components of G14.33— 0.63 are obtained to be 
U s = 11 ± 10 km s _1 toward the Galactic center at the 
source position, V s = — 1 ± 11 km s _1 in the local direction 
of the Galactic rotation, and W s = —4 ± 11 km s _1 verti- 
cally out of the Galactic plane toward the north Galactic 
pole. 

Here the uncertainties of 10 — 11 km s _1 in the derived 
peculiar motion are directly due to the uncertainties for 
the proper motion and radial velocity of G14.33— 0.64. 
The contribution from uncertainties in the Galactic con- 
stants i?o and 80 are negligible, because the Galactic 
rotation term is almost canceled out in the differential 
calculation. If we adopt the IAU standard values of 
Rq = 8.5 kpc and 80 = 220 km s™ 1 instead, the re- 
sulting peculiar motion becomes U s = 12 ± 10 km s _1 , 
y s = -l±llkms" 1 , and W s = -A± 11km s _1 . Therefore, 
the peculiar motion of G14.33— 0.64 is not significant in 



the direction of Galactic rotation (V s ) or in the direction 
out of the Galactic plane iW s ). For the source location of 
G14.33— 0.64 relative to the Sun in the Galaxy, the larger 
peculiar velocity component of G14.33— 0.64 toward the 
Galactic center (U s ) reflects a radial velocity larger than 
expected from the circular rotation model, which has led 
to the larger kinematic distances derived in the previous 
studies. Overall, G14.33— 0.64 shows no significant pecu- 
liar motion and is consistent with the circular Galactic 
rotation model. 
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Fig. 2. Parallax and absolute proper-motion measurements for G14.33— 0.64. Filled and open circles show positional evolution of 
mascr features 1 and 4, respectively, with respect to the reference positions at origin, (a) East offset (X) in mas from the reference 
positions of RA(J2000)=18 h 18 m 54. s 674440 for maser feature 1 and RA(J2000)=18 h 18 m 54. s 653410 for feature 4, as a function of 
time in days since the first epoch. Best-fitting models for parallax and proper motion are shown in solid curves and gray lines, 
respectively. Additional shifts are given for clarity: AX = +6, 3.5, 1, —1.5 mas for la, lb, lc, Id; —5, —7, —9 mas for 4a, 
4b, 4c, respectively, (b) North offset (Y) in mas from the reference position of Dec(J2000)= — 16°47'50" 26400 for feature 1 and 
Dec(J2000)= — 16° 47'50." 06500 for feature 4, as a function of time in days since the first epoch. Additional shifts are given for 
clarity: AY = +9, 6.5, 4, 1.5 mas for la, lb, lc, Id; —9, —11.5, —14 mas for 4a, 4b, 4c, respectively. 
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Table 3. Relative proper motions. 
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(1) Feature/spot ID. 'w-mcan' notates error-weighted mean of all spots of each feature and 'u-mcan' refers to 
unweighted- mean of all features (see text). (2) LSR velocity of the maser spot. (3) Number of detected epochs. (4) 
Detected epochs. (5) (6) Positional offset in mas toward east (X) and north (Y), respectively, from the reference 
spot 4b expected from the linear fit at the first epoch. (7) (8) Relative proper motion of the spot in X and Y, 
respectively, with respect to spot spot 4b. Estimated errors are shown in parentheses. 
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Fig. 4. (a) (b) Radio maps of G14.33— 0.64 showing H2O (dot) and methanol (triangle) maser positions superimposed on contours for 
6-cm continuum emission obtained from VLA archive data (program AH361). The angular resolution is 3" for the VLA observation 
(Hughes & MacLeod 1994) and the beam size is shown in gray at the left corner of each panel. Image noise level is <r = 0.07 mjy/bcam, 
and contours are linearly spaced and correspond to 4cr, 6cr, 8<r, ■ ■ ■, 22a. Peak intensity is 1.6 mjy/beam. Map origin is at the IRAS 
source position of 02000 =18 h 18 m 53. s 9, &2000 = — 16°47'39". Dots represent our VLBI absolute positions of H2O maser features in 
table 3. Numbers correspond to feature IDs in table 3. Our absolute position errors essentially come from errors of the reference 
quasar position J1825— 1718, which are 1.23 mas in RA and 1.97 mas in Dec (Fomalont et al. 2003). Triangles show the positions of 
44-GHz methanol masers mapped with the VLA by Slysh et al. (1999) with position errors of 0"2. Colors indicate the LSR velocity 
of the spots for both methanol and H2O maser emission, (c) Absolute proper motions of H2O maser features without correction of 
the Solar motion and Galactic rotation. Map origin (reference spot 4b) is at a 2 000 =18 h 18 m 54. s 653181 and 5 2 ooo = -16°47'50"07668 
(i.e., the position of spot 4b at the first epoch), (d) Internal motions of all maser features with the mean motion of the features of 
(fTx, /iV)=(0.95, —2.50) mas yr — 1 removed (without correction of the Solar motion and Galactic rotation). Map origin is the same 
as in figure 4c. A proper motion of 1.00 mas yr _1 corresponds to a linear velocity of 5.31 km s — 1 at a source distance of 1.12 kpc. 
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Fig. 5. Model of the Galaxy by Gcorgclin & Gcorgclin (1976), overlaid with the modified model by Taylor & Cordes (1993) shown 
in gray. '0' indicates the location of the Sun and 'GC the position of the Galactic center. The red square shows the new location 
of G14.33— 0.64 based on our parallax measurements, while the yellow square is the previously estimated position of G14.33— 0.64 
based on kinematic distances. Three star-forming regions, G35.20— 0.74 (blue diamond), G35.20— 1.74 (pink triangle) and W51 IRS2 
(green hexagon), possibly belonging to the Sagittarius spiral arm are also indicated with parallactic distances measured by Zhang 
et al. (2009) and by Xu et al. (2009) with the VLBA for 12-GHz methanol maser emission. Errors for all parallactic distances are 
also shown, which are within the size of dots for G14.33— 0.64 and for G35.20— 0.74. 




(deg) (j (deg) 

Fig. 6. Fits for the pitch angle of the Sagittarius spiral arm. The logarithm of Galactocentric radius R (measured in kpc) is plotted 
against Galactocentric longitude f} (in degrees). The Sun-center distance of 8.5 kpc was adopted, (a) G14.33— 0.64 (red square), 
G35.20— 0.74 (blue diamond), G35.20— 1.74 (pink triangle) and W51 IRS2 (green hexagon) arc plotted with parallaxes and associated 
uncertainties from this study, Zhang et al. (2009) and Xu et al. (2009). Gray lines show the best-fit straight lines from unweighted 
linear least-squares fitting to the data. The pitch angle i is obtained by taking the negative of the arctangent of the line slopes. 
(Note that we need to express ln_R in natural logarithm and /3 in radians to calculate the pitch angle.) Line A shows the fitting 
result from G14.33-0.64, G35.20-0.74 and G35.20-l.74, while line B is from G14.33-0.64, G35. 20-0.74 and W51 IRS2. (b) Same 
as (a), but with five sources (cyan dots) in the Local (Orion) arm (spur) also shown with precise parallax measurements. Line C is 
an unweighted straight line fit to W51 IRS2 and the five sources in the Local arm (see text). 



